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Being an atomically thin material, graphene is known to be extremely susceptible to its environment, 
including defects and phonons in the substrate on which it is placed as well as gas molecules that surround it. 
Thus, any device design using graphene has to take into consideration all surrounding components, and 
device performance needs to be evaluated in terms of environmental influence. However, no methods have 
been established to date to readily measure the density and distribution of external perturbations in a 
quantitative and non-destructive manner. Here, we present a rapid and non-contact method for visualizing 
the distribution of molecular adsorbates on graphene semi- quantitatively using terahertz time- domain 
spectroscopy and imaging. We found that the waveform of terahertz bursts emitted from graphene-coated 
InP sensitively changes with the type of atmospheric gas, laser irradiation time, and ultraviolet light 
illumination. The terahertz waveform change is explained through band structure modifications in the InP 
surface depletion layer due to the presence of localized electric dipoles induced by adsorbed oxygen. These 
results demonstrate that terahertz emission serves as a local probe for monitoring adsorption and 
desorption processes on graphene films and devices, suggesting a novel two-dimensional sensor for 
detecting local chemical reactions. 



While the highly tunable conductivity of charge carriers through gating is one of the attractive features of 
graphene^ it also means that the conductivity of graphene can be adversely affected by accidental 
doping in a rather uncontrolled fashion. It has been shown that adsorption of environmental gas 
molecules (even a single molecule) can drastically affect graphene's DC^ and opticaP conductivities. For detection 
and identification of molecular adsorbates on graphene, it is necessary to develop a versatile and convenient 
method that probes local conductivity or doping but does not require the fabrication of electrical contacts. Raman 
spectroscopy is one of the most useful tools for characterizing basic properties of graphene samples^, including the 
defect density, the degree of strain, and the number of layers. However, it does not provide direct information on 
the effect of molecular adsorbates, and their spatial distribution, on the conductivity. 

In the present study, we demonstrate visualization of the distribution of local molecular adsorbates on gra- 
phene using terahertz (THz) emission from a graphene-coated semiconductor surface. Specifically, we placed a 
large-area graphene film, grown by chemical vapor deposition^, on the (100) surface of InP, which is known to be 
an efficient source of coherent THz waves when excited by a femtosecond optical pulse^; see Fig. la for the 
experimental geometry. The waveform of the emitted THz burst sensitively changes due to the charge transfer 
between graphene and adsorbed molecules. According to classical electromagnetic theory, the electric field (E) of 
the generated THz wave is given by the time derivative of the transient current (/) induced by the optical pulse 
excitation, i.e., £thz ^ djidf. This equation indicates that the value and direction of the photoexcited current can 
be measured through the amplitude and polarity of the THz waveform, which are affected by the presence of 
molecular adsorbates. 

Results 

Figure lb shows how the THz waveform evolves in time under continuous excitation by femtosecond optical 
pulses from a Ti:sapphire laser, which remove gas molecules from graphene. Initially, there is a dip at —2.0 ps and 
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Figure 1 | Time-domain waveforms of terahertz radiation emitted from graphene-coated InP under different conditions, (a), A schematic diagram of 
the experimental geometry, (b), Time evolution of THz waveforms from graphene-coated InP under continuous excitation by femtosecond near-infrared 
pulses from a Tiisapphire laser. The curves are intentionally shifted vertically, (c), THz waveforms from InP (no graphene) under continuous excitation 
by femtosecond near-infrared pulses from a Ti:sapphire laser. The waveform remains the same even after ten minutes of continuous excitation, (d), THz 
waveforms from graphene-coated InP in vacuum under continuous excitation by femtosecond near-infrared pulses from a Ti:sapphire laser The 
waveform remains the same even after ten minutes of continuous excitation, (e), Time evolution of THz waveforms from graphene-coated InP under 
continuous excitation by femtosecond near-infrared pulses from a Tiisapphire laser. The vacuum pump is turned on at 10 minutes and the waveform 
changes drastically. 



a peak at —2.5 ps. After ~4 minutes of Ti:sapphire excitation, the 
2 ps dip vanishes, and a new positive peak appears at ~ 1.8 ps; at the 
same time, the amplitude of the 2.5 ps peak gradually decreases 
during 10 minutes of excitation, and the waveform does not change 
further. However, the waveform recovers its initial shape after the 
sample is exposed to air for a few hours without laser excitation, i.e., 
re-adsorption occurs. On the other hand, the waveform of THz 
radiation from a bare InP surface remains unchanged with time, as 
shown in Fig. Ic. Furthermore, the waveform of THz radiation from 
the graphene-coated InP (100) is also stable in time if the sample is 
placed in vacuum (—10"^ Torr), as demonstrated in Fig. Id. Finally, 
Fig. le shows dramatic time evolution of THz radiation from 
graphene-coated InP. The vacuum pump is turned on at 10 minutes, 
and the 2 ps feature suddenly changes its sign. All these results 
indicate that the change in the THz waveform with time arises from 
adsorption and desorption of atmospheric gas molecules on gra- 
phene induced by femtosecond laser excitation. 



In order to provide more insight into the origin of the THz wave- 
form changes due to molecular adsorption and desorption, we per- 
formed THz emission experiments under various conditions. 
Figures 2a, 2b, and 2c show the time evolution of the waveform of 
THz radiation from graphene-coated InP (100) in air, vacuum, and 
oxygen, respectively. Figure 2d shows the case where the sample was 
baked at 120° C under vacuum for 1 .5 hours to remove residual water 
at the interface between graphene and InP^'^ prior to exposing to 
oxygen. Before introducing gases, we kept the sample in a vacuum 
of 10"^ Torr. In all cases, after 10 minutes of THz emission measure- 
ment (Region I), we illuminated the sample by ultraviolet (UV) light 
(365 nm) for 5 minutes (Region II), motivated by the idea that UV 
illumination promotes oxidation of graphene^""^^. In Region I, the 
THz waveforms from all samples are unchanged in different envir- 
onmental gases and show similar characteristics, i.e., only the first 
peak can be seen. This result indicates that gas molecule adsorption 
dose not proceed much during 10-minute exposure to gases once 
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Figure 2 | Time evolution of THz emission from graphene- coated InP in 
different environmental gases, (a), Air, (b), vacuum, (c), oxygen, and (d), 
oxygen. In (d), the sample was pre-annealed to remove water before 
exposure to oxygen. In all cases, after 10 minutes of continuous THz 
emission measurements (Region I), the sample was illuminated by 
ultraviolet UV light (365 nm) for 5 minutes (Region II). After the UV light 
is switched off (Region III), the THz waveform slowly tends towards the 
original form. 

molecules are desorbed through vacuum pumping. Indeed, drastic 
changes are seen as a result of UV illumination, especially in air 
(Fig. 2a) and oxygen (Fig. 2c); the first peak at 2.0 ps suddenly 
becomes negative and the second peak at 2.5 ps suddenly appears 
upon UV illumination. On the other hand, little change is observed in 
vacuum (Fig. 2b). 

After the UV light is switched off (Region III), the THz waveform 
in air and oxygen slowly returns towards the original shape but does 
not completely recover the initial waveform in the time frame of the 
measurements. Basically no similar change is observed when the UV 
light is turned off in the vacuum case (Fig. 2b). These results suggest 
that the existence of oxygen molecules is necessary for changes to 
occur in the waveform of THz radiation emitted from graphene- 
coated InP (100) upon UV illumination. Finally, for the pre-annealed 
sample measured in the presence of oxygen (Fig. 2d), the THz wave- 
form has one peak in Region I and shows a smaller change under UV 
light illumination than the case without pre-annealing (Fig. 2c). In 
Region III in Fig. 2d, the waveform completely recovers the initial 
waveform whereas in Fig. 2c the waveform does not recover the 



original. These results show that the dramatic UV response is due 
to the co-existence of oxygen and water adsorbates on graphene, 
consistent with previous results reported by Mitoma et aU^. 
Moreover, adsorbed oxygen molecules, after UV illumination in 
the presence of water, are not easily removable by femtosecond laser 
excitation (Region III in Figs. 2a and 2c). 

Discussion 

The emission mechanism of THz radiation from an InP surface is 
primarily the current surge effect due to the built-in electric field in 
the surface depletion layer^. We believe that the observed changes in 
THz waveforms are caused by the existence of electric dipoles 
induced by adsorbed oxygen molecules, which modify the surface 
electric field. The schematic band diagrams in Fig. 3 show the 
expected flow of photoexcited carriers both in the (a) absence and 
(b) presence of adsorbed oxygen. The Fermi level of bulk semi- 
insulating InP (100) is located 0.75 eV below the conduction band 
bottom while at the surface it is pinned at 0.45 eV below the con- 
duction band bottom^'^^. Thus, the band bends downward toward the 
surface, and the photoexcited current driven by the surface field flows 
toward the inside of the InP (100) substrate (Fig. 3a). Judging from 
the fact that the THz waveform is identical between bare InP (100) 
and graphene-coated InP (100) under vacuum, it can be assumed 
that pristine graphene alone, without adsorbates, does not change the 
surface band bending. On the other hand, exposure to air for an 
extended period results in a charge transfer between graphene and 
adsorbates^'^^"^^. Theoretical calculations suggest that, while gra- 
phene is hydrophobic, it is still capable of adsorbing small amounts 
of oxygen and water molecules^^"^^. Moreover, it has been reported 
that the shift of the Fermi level induced by gas adsorption strongly 
depends on the substrate^^'^^'^^. Therefore, it is likely that adsorbed 
gas molecules predominantly exist between graphene and InP (100), 
as shown in Fig. 3b, forming local electric dipoles at the InP surface, 
which induce changes in the waveform of THz radiation. Specifically, 
the electric field generated by the local dipole between graphene and 
InP increases the surface potential of InP (100), which tends to bend 
the band upward toward the surface. As a result, the photoexcited 
current flows toward the surface (Fig. 3b), opposite to the case with- 
out molecules (Fig. 3a). 

In general, the waveform of THz radiation emitted from 
graphene-coated InP should be a superposition of two waveforms 
representing, respectively, the band diagrams in Figs. 3a and 3b. The 
relative weights of the two components should depend on the density 
of local dipoles within the focal spot of the excitation laser beam, 
which gradually decreases with time under laser excitation. We 
explain our theoretical model in a semi- quantitative manner with 
experimental data in Fig. 4. Figure 4a shows two different waveforms 
- the one taken under high vacuum (red line, representing Fig. 3a) 
and the other taken after exposure to air for a few days (black line, 
representing Fig. 3b). The observed THz waveforms (Fig. 4c) for 
graphene-coated InP in air at different times are simulated with 
waveforms consisting of a superposition of these two observed wave- 
forms in Fig. 4a at different relative weights. Figure 4b shows the 
relative weights of the waveforms taken after exposure to air used for 
the simulation of experimental data. The simulated waveforms 
(Fig. 4d) agree with the observed waveforms very well, with appro- 
priately changing relative weights in the superposition. Furthermore, 
the UV-induced changes in THz waveforms can also be explained by 
our model. UV light promotes the oxidation of graphene in the 
presence of oxygen and water molecules^^"^^; this reaction is rapid 
and creates tightly-bonded chemisorbed oxygen. Therefore, desorp- 
tion of oxygen molecules via femtosecond laser excitation becomes 
more difficult, and the THz waveform does not fully recover the 
initial one in Region III (Figs. 2a and 2c). 

Finally, Fig. 5a shows a map of the amplitude of THz radiation 
emitted from graphene-coated InP (100) after exposure to air for a 
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Figure 3 | Band diagrams explaining the sign change of THz amplitude due to adsorbed oxygen molecules, (a), Band diagram for pristine graphene 
on InP. The photo excited current flows towards the substrate, (b), Band diagram for graphene on InP with adsorbed oxygen molecules between them. 
The charge transfer between graphene and oxygen creates dipoles, which modif)^ the band bending. The photoexicted current now flows towards the 
surface. 



few days. The amplitude of the first peak greatly varies, exhibiting 
both positive and negative values; a reference InP substrate without 
graphene shows almost constant values (see Supplementary Infor- 
mation). The THz waveforms corresponding to Point 1 and Point 2 
in the map (Fig. 5a) as wefl as that for a bare InP (100) substrate are 
shown in Fig. 5b. The waveform at Point 1 is simflar to that from a 
highly polarized surface due to adsorption of oxygen molecules, and 
the waveform at Point 2 is qualitatively the same as that from the bare 
InP. These results indicate that this method using THz emission is a 



valuable tool for measuring the distribution of local dipoles induced 
by oxygen adsorption. 

Methods 

Terahertz time-domain spectroscopy. The graphene-coated InP (100) and the 
reference bare InP crystals were optically excited by near-infrared pulses with a 
duration of — 100 fs and photon energy of 1.55 eV generated by a Ti:sapphire laser. 
The excitation laser beam was delivered to the sample surface at an incident angle of 
45°. The power of the excitation laser beam was set at 20 mW. The diameter of 
excitation laser spot was —100 |im, and the fluence was three orders of magnitude 
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Figure 4 | Laser- induced desorption dynamics of oxygen molecules from graphene observed through THz waveforms, (a), Two kinds of waveforms 
emitted from graphene-coated InP, one taken under high vacuum (red line) and the other taken after exposure to air for several days (black line), (b), The 
time dependent relative weights of the waveforms taken after exposure to air used for the simulation of experimental data, (c), THz waveforms for 
graphene-coated InP in air at different times, (d), Simulated waveforms consisting of a superposition of two observed waveforms of different 
relative weights. The simulated waveforms agree with the observed waveforms very well, with appropriately changing relative weights in the 
superposition. 
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Figure 5 | THz imaging of adsorbed oxygen on graphene. (a), Amplitude 
mapping of THz radiation emitted from graphene-coated InP after 
exposure to air for a few weeks, (b), THz waveforms corresponding to 
Point 1 and Point 2 indicated in the map as well as that for a bare InP 
substrate. The waveform at Point 1 indicates a highly polarized surface due 
to adsorption of oxygen molecules, and the waveform at Point 2 is 
qualitatively the same as that from the bare InP. 

lower than the damage threshold previously reported^^. THz radiation was mainly 
emitted from the surface of the InP (100) substrate through ultrafast surface surge- 
currents with the range of the fluence in this experimet^'^*^. THz radiation was emitted 
into free space and then focused onto a dipole-shaped low-temperature-grown GaAs 
photoconductive switch via a pair of off-axis parabolic mirrors. The output of the 
photoconductive antenna was fed into a lock-in amplifier, referenced to a 2.0 kHz 
optical chopper signal. THz radiation images were obtained by raster- scanning of an 
x-y mechanical stage on which the sample was mounted. For UV light illumination, a 
light emitting diode with 365 nm wavelength was utilized. The spot diameter was 
~8 mm, and the power was 36 mW. The gases used in the THz emission 
measurements were ambient air, dry air, 99.9995% pure N2, and 99.999% pure O2. All 
measurements in atmospheric gases were performed at 1 atm. 

Sample fabrication. The studied graphene samples were synthesized on copper (Cu) 
foil using hexane as a liquid precursor in a chemical vapor deposition system^. After 
growth, a thin poly-methyl-meth-acrylate (PMMA) film was deposited on the 
graphene/Cu substrate as supporting layer for the graphene transfer. The underlying 
Cu substrate was then dissolved in dilute HNO3 and washed several times in 
deionized water. The floating graphene/PMMA film was transferred onto an InP 
substrate. The sample was immersed into acetone solvent to remove PMMA, and then 
dried in critical point dryer to avoid damage of the partially deposited graphene layer 
by the local capillary forces arising in conventional drying processes. 
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